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Abstract. We present Doppler Images of the RS CVn binary system HR 1099 (V711 Tau) from spectra taken in 
two different sites, KPNO and MSO, during the MUSICOS 1998 campaign. Contemporaneous APT photometry 
is used to constrain the Doppler Images. The resulting maximum entropy reconstructions based on the least- 
squares deconvolved profiles, derived from ~2000 photospheric absorption lines, reveal the presence of starspots 
at medium-high latitudes. We have obtained maps for both components of the binary system for the first time. 
The predominant structure in the primary component is an off-centered polar spot, confirming previous works 
on the same target by using independent codes. The result is verified by using both data sets independently. The 
lower spectral resolution data set gives a less detailed map for the MSO data set. The images obtained for the 
secondary component show a low latitude spot around orbital phase 0.7. This spot seems to mirror the structure 
seen on the primary. It might suggest that tidal forces may influence the spot distribution on this binary system. 
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1. Introduction 

RS CVn binary systems consist of a chromospherically 
active evolved star tidally locked to a main-sequence or 
sub-giant companion. Short orbital periods ranging from 
a few days to 20 days are typically observed. The RS 
CVn high level of activity has been measured from ra- 
dio to hard X-rays and, for fast rotators, it approaches 
the saturation limits for chromospheric, transition region 
and coronal emission. One of the striki ng aspects of these 
systems is their propensi ty to flare l|Dovle et alJ |l992; 
iGartia- Alvarez et alJ l2003h Moreover, RS CVns show op- 
tical photometric variations which are believed to arise 
from the rotational modulation of ph otospheric spots 
ijRodono et alJ Il98ri iBvrne et al.lll99oj) . large scale ver- 
sions of dark solar spo ts, that p rovide evide nce for large- 
scale magnetic fields. iDonatil l)l999j) and iDonati et alJ 
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(2003) showed that Zeeman Doppler Imaging can also pro- 
vide evidence of such magnetic structures. Short-period 
RS CVn-like systems, through their rotational modula- 
tion, can then provide information on the morphology and 
two-dimensional structure of active regions in stellar at- 
mospheres. 

iDeutscbl (Il958lll97(t) first proposed the use of line pro- 
files a s a tool to map stellar surfaces. iGoncharskv et all 
( 1977) developed the fi rst inversion technique with reg- 
ularization, but it was IVoet k. Penrodl {l983) who first 
called it Doppler Imaging. The technique exploits the cor- 
respondence between the wavelength position across a ro- 
tationally broadened spectral line and the spatial posi- 
tion a cross the stellar disk at a given rotation phase. Vogt 
(1987) addressed the difficulties in inverting the spectral 
line profiles into an image of the stellar disk. Doppler 
imaging has provided one of the few direct pieces of ev- 
idence for the existence of starspots, or localised plages, 
but the technique is severely restricted. The data must be 
of very high quality (S/N > 200) and the method can only 
be applied to rapidly rotating stars (wsini > 20 kms -1 ) 
up to a limit where the signal is swamped by terrestrial 
and stellar line blends (usini > 160 kms^ 1 ). The stel- 
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lar inclination is limited to a range of 20°<i<90°, and 
the physical properties of the atmosphere must also be 
known, although the process is fairly insensitive to the de- 
tailed physics of spectral line formation. Several reviews, 
explaining the principl es of this technique, can be found 
in the literature (e.g., iRice et al1ll989t ICollier Cameronl 
l200lllRicdl2r)03) . 

The target under consideration, HR 1099 (V711 Tau, 
03 h 36 m 47 s +00°35'16" (J2000), V=5.64, B-V=0.92), is 
a close double-lined spectroscopic binary. At a dis tance 
of 29 pc (The Hipparcos and Tycho Catalogues (ESA 
I1997IU . HR 1099 is the nearest and brightest of the clas- 
sical RS CVns with a K1IV primary and a G5V sec- 
ondary tidally locked in a 2"? 8 orbit. iFekell l|l983h pro- 
vided the orbital parameters of the binary and obtained 
masses, radii, and spectral types for both components. 
IGarcia- Alvarez et alJ l|2003|) obtained the most recent or- 
bital parameters of the system. The K sub-giant nearly 
fills its Roche lob e and is by far the more active of the 
two components ( Avres fc Linskvl Il982t iRobinson et alJ 



presence of a long-lived (>11 yr) polar spot, together 
with transient (<1 yr) lo w-latitude spots on the sur- 
face of the active K star. IStrassmeier fe Bartusl l|200fjft 
reported latitudinal and longitudinal migration of indi- 
vidual starspots and their morphological evolution. They 
confirme d the poleward migrati on scenario previously de- 
tected bv lVoet fe HatzeJ (1996). We note that all the au- 
thors, except Gondoin (1986), have found a polar spot on 
the surface of HR 1099. Strong magnetic fields of the order 
of 1000 G have also been detected on HR 1099 using; the 
Zeeman-Doppler I maging technique I Donati et al.1 Il990[ 



Evolutionary models suggest that mass transfer 
from the K primary onto the G secondary may be- 
gin within 10 7 years IjFekell 1 198311 . The primary exhibits 
conspicuous signatures of chromospheric activity, such 



as strong and variable Call H & K and Ha emission 


jRodono et al. 1987llDemDsev et s 


1. 1996: Robinson et alJ 


1996: iGarci'a-Alvarez et al.l 


2003 


. This is due to its 



tidally induced rapid rotation, combined with the deep- 
ened convec t ion zo ne of a post-main-sequence envelope. 
lAvres et alJ l|200ll) reported that the Kl IV component 
dominates the X-ray emission in the binary system HR 
1099. 

The G5 dwarf is probably also active, having a suf- 
ficiently deep convection zone and fast rotation to host 
an efficient dynamo. In particular, due to the presence 
of spots, main-sequence G-type stars with rotation peri- 
ods of about 3 days are expected to show V-band light 
curve amplitudes up to 0.10 mag, as can be inferred 
fro m the rotation-activi ty relations at photospheric lev- 
els ijMessina et al.ll200l[l . However, because of the small 
luminosity ratio in the V-band (Lqsv/Lkiiv — 0.27), the 
contribution from the G5 component to the observed mag- 
netic activity manifestations in the optical band is signif- 
icantly smaller than from the Kl IV primary. 

Extensive broadband photometric monitoring has been 
a key factor in verifying the presence of evolving spots 
on the surface of HR 1099 and deducing basic spot pa- 
rameters (i.e., size, temperature, longitude). HR 1099 is 
also among the brightest spotted stars and, therefore, 
a imp ortant target for Doppler Imaging. IVogt fc Penrodl 
ll 19831) p r esente d the first Doppler Image of HR 1099. 
iGondoir] l(l986h used a technique, originally applied to 
equivalent width (hereafter EW) variations of Ap stars, 
to map the HR 1099 surface. iFoing et all <|l994^ obtained 
complete phase coverage for Doppler imaging o f spots dur- 
ing th e MUSICOS 1989 campaign on HR 1099. IVogt et alJ 
l)l999ll presented a seminal paper including 23 Doppler 
images obtained from 1981 to 1992. They showed the 
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Th is work follows that done by IGarcia- Alvarez et alJ 
l|2003h on HR 1099 and based on the MUSICOS 1998 cam- 
paign. Using photometry spot modelling, they obtained 
maps of the distribution of the spotted regions on the pho- 
tosphere of the binary components. Those maps were de- 
rived using the Maximum Entropy method and Tikhonov 
regularization criterion on photometric data sets only. The 
maps show the Kl IV primary to be the more active with a 
large spotted region in the northern hemisphere, centered 
around phase 0.85 (computed according to the cphemeris 
given in Eq. (1)). The detection of optical and X-rays 
flare events around the same phase is in good agreement 
with a spa tial link between flare s and a ctive regions as re- 
ported bv IGarcia- Alvarez et alJ l|2003ri . In this paper, we 
present the results of the Doppler Imaging technique ap- 
plied to the spectroscopic data sets obtained at Kitt Peak 
National Observatory (hereafter KPNO) and Mt. Stromlo 
Observatory (hereafter MSO), during the MUSICOS 1998 
campaign on HR 1099. Observations and data analysis are 
described in Sect. 2; in Sect. 3 we present the image recon- 
struction process, and in Sect. 4 the results and discussion 
are given. 

2. Observations and Datasets 

The analysis of the rotational modulation of the opti- 
cal bands flux (i.e., U, B, V) allows us to obtain maps 
of the stellar photosphere or more precisely, tempera- 
ture (in the case of multiband photometry), total area 
and longitudinal distribution of the spotted regions. The 
deeper the minimum in the light curve, the larger the 
spotted region on the stellar disk for a given temper- 
ature difference between the unperturbed photosphere 
and the spots. The shape of the light curve is usually 
not sinusoidal, indicating that more than one active re- 
gion i s required to interpr et the optical flux modulation 
(e.g., iRodono et alJ fl9861 . However, apart from eclips- 
ing binaries, no reliable constraints can be derived from 
photomery for what concerns the latitudinal distribution 
of spotted regions. Moreover, due to limb-darkening ef- 
fects only spots located at low latitudes in equator-on 
stars, or more generally, near the substellar latitudes con- 
tribute appreciably to the observed light curve modula- 
tion. Thus, while photometry is most sensitive to spots at 
latitudes much higher/lower than the substellar ones, the 
spectroscopy reconstructions have the best latitude dis- 
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Table 1. The log of observations of HR 1099 and Standards 
from KPNO and MSO during the MUSICOS 1998 campaign. 



KPNO 



Object 


UT Start 


Exp. time 


No. of 


Comments 




1998 Nov 24 


(sec.) 


frames 




HD 218045 


01:16 


300 


2 


B9III telluric std. 


HD 5268 


01:33 


300 


1 


G5 IV template 


HR 1099 


01:53 


1200 


4 


Target Star 


Object 


UT Start 


Exp. time 


No. of 


Comments 




1998 Nov 25 


(sec.) 


frames 




HD 218045 


01:23 


300 


1 


B9III telluric std. 


HD 5268 


01:32 


300 


1 


G5 IV template 


HR 1099 


01:51 


1200 


4 


Target Star 


HD 82074 


11:44 


300 


1 


G6 IV template 


HD 92588 


12:01 


300 


1 


Kl IV template 


Object 


UT Start 


Exp. time 


No. of 


Comments 




1998 Nov 26 


(sec.) 


frames 




HD 218045 


05:23 


300 


2 


B9III telluric std. 


HR 1099 


06:03 


1200 


5 


Target Star 


Object 


UT Start 


Exp. time 


No. of 


Comments 




1998 Nov 27 


(sec.) 


frames 




HD 218045 


02:29 


300 


2 


B9III telluric std. 


HR 1099 


02:56 


1200 


7 


Target Star 


MSO 


Object 


UT Start 


Exp. time 


No. of 


Comments 




1998 Nov 29 


(sec.) 


frames 




HR 1099 


12:16 


(100 


3 


Target Star 


HD 36079 


16:11 


300 


1 


G5 II template 


Object 


UT Start 


Exp. time 


No. of 


Comments 




1998 Nov 30 


(sec.) 


frames 




HD 5112 


09:48 


300 


1 


M0 III template 


HR 1099 


10:10 


000 


4 


Target Star 


HD 39523 


17:51 


300 


1 


Kl III template 


Object 


UT Start 


Exp. time 


No. of 


Comments 




1998 Dec 01 


(sec.) 


frames 




HD 218045 


09:43 


300 


1 


B9 III template 


HD 886 


09:51 


300 


1 


B2IV telluric std. 


HR 1099 


10:14 


800 


4 


Target Star 


HD 39523 


17:57 


300 


1 


Kl III template 


Object 


UT Start 


Exp. time 


No. of 


Comments 




1998 Dec 02 


(sec.) 


frames 




HD 218045 


09:44 


300 


1 


B9III template 


HR 1099 


10:17 


000 


1 


Target Star 


Object 


UT Start 


Exp. time 


No. of 


Comments 




1998 Dec 03 


(sec.) 


frames 




HD 5112 


10:55 


300 


1 


M0 III template 


HR 1099 


11:16 


800 


5 


Target Star 


HD 36079 


16:53 


300 


1 


G5 II template 



trimination at higher latitudes t han substellar ones (e.g., 
lUnruh fc Collier Cameror]ll997l) . Therefore, in order to 
obtain the best possible reconstruction of the distribution 
of brightness on a stellar photosphere it is necessary to 
combine both spectroscopic and photometric data sets. 



2.1. MUSICOS 98 campaign: Optical Spectroscopy 

The data analysed in this work were obtained at two sites, 
namely: Kitt Peak National Observatory (KPNO), USA 
and Mt. Stromlo Observatory (MSO), Australia, as part 
of the MUSICOS 1998 campaign. 

The KPNO data set was obtained on 1998 November 
24-27 using the coude spectrograph and the 3kxlk F3KB 
CCD at the 0.9m coude feed telescope. The wavelength 
range of the extracted spectra is 5465-7060 A, recorded 
over 23 orders with a resolving power of R=65 000. A to- 
tal of 20 spectra of exposure time 1200 s were obtained, 
along with spectral standard stars with effective temper- 
atures similar to those of the HR 1099 K1IV and G5V 
components. 

The MSO data set was obtained on 1998 November 29 
- December 3 using the echelle spectrograph on the 74-inch 
telescope. The 4kx2k CCD17 chip was used. A bluer set- 
ting with a wavelength coverage of 4540-6700 A, recorded 
over 43 orders with a resolving power of R=35 000, was 
used at this site. A total of 17 spectra of exposure time of 
800 s/900s were obtained, together with spectral standard 
stars with effective temperatures similar to the HR 1099 
components. 

The details of the observations are recorded in Table 1. 
A more detailed log of observations of the target HR 1099, 
during the MUSICOS 1998 campaign, can be found in 
iGarcfa- Alvarez et al 1 ll2003h . 

The spectra were extracted using the standard reduc- 
tion procedures in the IRAF 1 packages. Background sub- 
traction and flat-field correction using exposures of a tung- 
sten lamp were applied. The wavelength calibration was 
obtained by taking spectra of a Th-Ar lamp. A first-order 
spline cubic fit to approximately 45 lines achieved a nom- 
inal wavelength calibration accuracy which ranged from 
0.06 to 0.11 A. The spectra were normalised by a low- 
order polynomial fit to the observed continuum. Finally, 
for the spectra affected by water lines, a telluric correction 
was applied. 

To produce a p hase-folded light curve we ad opted the 
ephemeris given bv lGarcia- Alvarez et alJ l)2003|) : 

HJD = 2451142.943+ 2. 83774S (1) 

2.2. Photometry 

During the MUSICOS 98 campaign, contemporary 
ground-based photoelectric observations were obtained 
with two different Automatic Photoelectric Telescopes 
(APTs): a) the 0.80-m APT-80 of Catania Astrophysical 
Obse rvatory (CAO) on Mt. Etna, Italy l|Rodono et alJ 
1200 lj) ; b) the 0.25-m Tl Phoenix APT of Fairborn 
Observatory at Washington Camp, AZ, USA, that is 

1 IRAF is distributed by the National Optical 
Astronomy Observatories, which is operated by the 
Association of Universities for Research in Astronomy, 
Inc., under cooperative agreement with the National 
Science Foundation. 
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Fig. 1. KPNO (left panel) and MSO (right panel) grey-scale time-series images of the deconvolved spectra of both components 
of HR 1099, for each site. The plots are residuals after subtraction of the mean profile for each data set (white features are spot 
transients). The vertical lines are at the vsini of the K1IV component. 



Table 2. Least-squares multiplex gain for each data set. The servations. This section describes the steps that have been 
output S/N ratios are calculated from the LSDprocess. followed to reach this aim. 



Data Set Input LSD output Gain Effective 



S/N ratio S/N ratio S/N ratio No. of lines 



KNPO 24th 


65 


1790 


27.7 


1711 


KNPO 25th 


73 


2038 


27.9 


1711 


KNPO 26th 


74 


2052 


27.8 


1711 


KNPO 27th 


78 


2163 


27.9 


1711 


MSO 29th 


86 


3656 


42.3 


2249 


MSO 30th 


92 


3861 


42.2 


2249 


MSO 01th 


83 


3539 


42.4 


2249 


MSO 02th 


59 


2489 


41.9 


2249 


MSO 03th 


87 


3634 


41.6 


2249 



managed a s a multiuser telescope l)Bovd et all Il984t 
ISeedslll99ffl . Both telescopes are equipped with stan- 
dard Johnson UBV filters and uncooled photoelectric 
photometers. Although repeated observations were sched- 
uled on each night during the campaign to obtain com- 
plete rotation phase coverage, bad weather conditions re- 
sulted in light curve gaps around phases 0.30 and 0.70. 
Details on the observation a nd reduction can be found in 
iGarcia- Alvarez et all l|2003|) . 

3. Image Reconstruction 

In order to obtain the best possible results, the Doppler 
imaging technique requires careful preparation of the ob- 



3.1. Continuum fitting 

The best possible fit to the observed continuum is re- 
quired, in order to obtain the best results from the least- 
squares deconvolution. The continuum of an echelle spec- 
trum can be estimated by using several methods. The 
blaze function dominates the shape of the continuum in 
the echelle orders. For CCDs of sizes of 2kx2k pixels, or 
even larger, the continuum is generally well approximated 
by a carefully chosen spline fit. The number of knot points 
required in the spline fit varies from one data set to the 
next. In this case, eight knot points were used for standard 
stars on both KPNO and MSO data sets, while eleven knot 
points were used for the target star, HR 1099. 

We have taken as a master frame a spectrum in which 
one component of the system is in front of the other. This 
allows us to observe the largest amount of continuum pos- 
sible. This master continuum was then divided by the fit to 
the continuum of the two slowly rotating standards that 
match the spectral types of the components of the HR 
1099 system. Note that different standards with different 
spectral types were used for the KPNO and MSO data 
sets, respectively. This contin uum fitting process is sim- 
ilar to that first described bv ICollier-Cameron fc Unruhl 
l|l994j) . This approach has the advantage that the lines 
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Fig. 2. Time-series profiles and maximum-entropy-regularized fits to the HR 1099 binary system data observed at KPNO (left 
panel) and MSO (right panel) during the MUSICOS 1998 campaign. The number to the right of each profile is the phase of the 
observation. We have plotted the spectra ordered in phase for easier comparison between them. The level of the continuum has 
been systematically shifted to simplify the comparison of all the profiles. 



do not suffer from rotational broadening, giving a greater 
number of continuum windows. Since the continuum level 
and shape of HR 1099 differ slightly from those of the pho- 
tospheric standards, a low order (quartic) polynomial was 
fitted to this divided spectrum to remove residual misfits, 
and to normalise the spectrum, resulting in a standard or 
master continuum fit. 

Normalisation of the HR 1099 spectra were carried out 
by dividing each target spectrum by the master continuum 



frame, in order to obtain an essentially flat residual. This 
residual was then fitted with a quartic polynomial to re- 
move residual misfits and any remaining tilt, mainly due to 
extinction effects as the star rises and sets. Slight changes 
in continuum tilt due to the starspot contribution are also 
removed in the process. This fit, when multiplied by the 
master continuum fit, results in the correct continuum for 
the target spectrum. 
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3.2. The local specific intensity profiles. 

A spotted star is considered to have two temperature com- 
ponents (i.e. photospheric temperature and spot temper- 
ature). In the case of a binary system, two photospheric 
temperatures and two spot temperatures, for the pri- 
mary and secondary component respectively, are needed. 
Two template stellar spectra to mimic the photosphere 
temperature were used, namely: a K1IV (HD 92588 and 
HD 39523 for the KPNO and MSO data set, respectively) 
with a photospheric temperature T=4800K for the pri- 
mary component, and a G5 V (HD 5268 and HD 36079 for 
the KPNO and MSO data set, respectively) with temper- 
ature T=5400K for the secondary component. The tem- 
plate stellar spectra MO III (HD 5112), with T=3 800K, 
was used for spot temperature in both components, in 
the MSO data set. The KPNO data set does not con- 
tain M type standard, so to mimic the spot temperature 
we used an artificial echelle spectrum generated by us- 
ing a r=3 800K li ne list (Vienna Atomic Line Database; 
iKunka et al1ll999l) . The choice of starspots with temper- 
ature T=3 800K is based upon similar single stars or 
RS CVn subgiants for which the spot umbrae is 500- 
1500K cooler than the photosphere. These values were 
obtained either by u sing multicolour photometry (i.e., 
lOlah et al.ll2002lll997l) or by using Ti O absorption bands 
H.e.. lSaar et a,lJl200lUNeff et al.lll995h . We used the same 
spot temperature in both binary components in order to 
simplify the analysis. However, if the spot temperature is 
wrong, incorrect limb-darkening coefficients will be used 
but, because of the small effective intensity of the spot 
(typically less than l/10ththe photospheric contribution), 
this will only affect the spot filling factor obtained in the 
images rather than the shape of any spot features. There 
are indications that the spot properties are sensitive to 
surface gravity, therefore, when the MO III spectra is used 
as template, for spots on the components of HR 1099, we 
should take into account that we are working with a cer- 
tain level of approximation. 

M dwarf spectra contain molecular absorption bands, 
whose opacities are poorly re presented in most available 
model atmospheres. However, iBarnes fc Collier Cameron! 
(2001) concluded that the use of an incorrect line list does 
not severely affect the shape of the deconvolved profile. All 
the template stellar spectra were deconvolved, using the 
method of least-squares deconvolution, in the same man- 
ner and using the same scaled continuum frame, ensuring 
correct equivalent width. 



3.3. Limb-darkening coefficient 

The centroidal wavelength was 6043 A and 5423 A for 
KNPO and MSO data sets, respectively. It was calcu- 
lated from lines in the line list, taking into account 
both the strength and the variance from the standard fit 
used in the spectral deconvolution procedure (hereafter 
SPDECON). Due to the fact that the limb-darkening co- 
efficient is a fairly linear function of wavelength, we used 



linear interpolation to obtain a valu e from local thermo- 
dynamic equilibrium (PT E) models ijDiaz-Cordoves et alJ 
Il995l IClaret et alJll995|) . The limb-darkening coefficients 
for the primary component K1IV are 0.761 and 0.729 for 
the 3 500 K spot temperature and for the 4 750 K photo- 
spheric temperature respectively. For the secondary com- 
ponent G5 V the limb-darkening coefficients are 0.710 and 
0.646 for the 3 600 K spot temperature and for the 5 500 K 
photospheric temperature respectively. In this analysis 10 
different limb angles (i.e., /i=0.1,0.2,...,1.0) for each tem- 
perature were used. 



3.4. Least-squares deconvolution 

We use the technique of least-squ ares deconvol u tion ( LSD 
hereafter), first implemented bv iDonati et alJ l|l997h . to 
derive an average line profile. This technique can be 
used to combine the information content of many weak 
and intermediate-strength photospheric absorption lines, 
by computing average profiles from thousands of spec- 
tral echellogram lines simultaneously. We assume that an 
echelle spectrum contains the same starspot signatures in 
all the photospheric Doppler-broadened profiles. Although 
the amplitude of the signature may change for different 
lines, the morphology should remain the same. A detailed 
mathem atical explanation o f th e LSD technique can b e 
found in lDonati etaD l|l997h and lCollier Cameron! l|200lh . 



The maximum achievable gain using LSD is approx- 
imately (number of lines) 1 / 2 , however this is an upper 
limit not reachable in most of the cases. In our KPNO 
and MSO data sets ~1711 and ~2249 photospheric ab- 
sorption lines are available, respectively. Table 2 shows 
the S/N obtained by our deconvolution code SPDECON 
along with other statistics. The velocity bin width in the 
deconvolved profiles is generally set to match the mean 
velocity bin size of each pixel on the CCD. In the case of 
the KPNO data set the deconvolved profiles for all the HR 
1099 observations were set to a bin width of 1.3 kms -1 , 
while for MSO the value was set to 3.0kms _1 . 

Lines in the regions of known strong TiO bands (such 
as at 6782, 7055 and 7126 A) were removed from the line 
list. Lines with wavelength >7200 A are not included in 
the deconvolution due to the presence of numerous telluric 
lines. All the regions around the strongest atomic lines 
such as Ha (6555-6570 A), H/3 (4840-4880 A), Na Doublet 
(5880-5910 A) and other lines that can be affected during 
flares such as He I D 3 (5870-5880 A) and Mg I (5160-5190 
A), were also removed from the line list. 

The time-series spectra obtained after applying least- 
squares deconvolution are shown in Fig. ^ I n the time- 
series spectra it is possible to identify a very clear starspot 
signature, visible in both data sets. This feature seems to 
reaches its maximum around phases 0.6 and 0.7, which 
is close to the phase of flux minimum in the photometric 
data set. We also note that, between phase 0.1 and 0.2, 
there is an apparent lack of starspots, which is in agree- 
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Fig. 3. Photometric light curve of HR 1099 obtained contemporary to the MUSICOS 98 campaign. Top panel: Filled dots 
represent differential V-band magnitudes and open dots the maximum entropy regularized fit. Bottom panel: Residuals. 

Table 3. Estimated system parameters from maximum-entropy minimized \ 2 fits to data from each site. 7 denote the systemic 
velocity, K\ is the velocity amplitude for the primary component, EW pr i and EW SC c are the equivalent width of the primary 
and secondary components respectively, MijMx is the mass ratio, i is the inclination angle and e is the eccentricity. 
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EW pri 


EWsec 


Period 


Rpri 


R-scc 


M2/M1 


i 
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(kms" 1 ) 


(kms -1 ) 


(A) 


(A) 


(d) 


(Ho) 


(Ro) 




(deg) 




KPNO 


-12.74 


49.66 


117 


139 


2.8377 


3.3 


1.1 


0.789 


40 





MSO 


-15.69 


49.66 


82 


210 


2.8377 


3.3 


1.1 


0.789 


40 






merit with the maximum in the photometric light curve 
around those phases. 



3.5. Determination of system parameters 

Accurate system parameters, such as the systemic veloc- 
ity, equivalent width, inclination angle, etc, are needed 
to obtain the best results from Doppler Imaging. In fact, 
spurious features are pr oduced as a result of incorrect in - 
put parameters values ijCollier-Cameron fc Unrublll994jh 
System parameters can be determined empirically by min- 
imization of either spot area (i.e., in the reconstructed 
images) or minimization of the x 2 fit to the data. We 
have chosen the \ 2 minimization method, since this is 
more satisfactory for spectra of high S /N ratios and does 
not introduce prob lems as the spot area method does 
(Bar nes et al.l200o|) . The \ 2 method works as follow: After 
obtaining initial estimates for the input parameters, we 
keep all them fixed, except one. A large number of max- 



imum entropy regularised iterations are then performed 
by our Dopplcr-Imaging code DoTS , in order to converge 
to a y 2 minimum. Th e initia l parameters were taken from 
iGarcia- Alvarez et al.l l)2003h . Table 3 shows the derived 
system parameters. There is good agreement between the 
values derived from the KPNO data set, and those inde- 
pendently derived from the MSO data set. Only the value 
of the axial inclination is poorly known. There have been 
several atte mpts to me asure it, yielding estim ates ranging 
from -33° (|Fekellll983f) to ~45° l|Donatll999|) . We choose 
i=40° since it gives the minimum chi-squares. Errors of 
±10° generally have an effec t only on the spot filling fac- 
tor as a function of latitude ^Barnes! 1999|) . However, very 
strong polar features could appear by overestimating the 
axial inclination, under-representing low-latitude features 
as well. The opposite is true for underestimation of the 
axial inclination. 

Note that the internal shift of the instruments, which 
could affect the radial velocity measurement, are corrected 
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Fig. 4. Maximum-entropy-regularized image reconstructions, for the primary component (Kl IV) of the HR 1099 system, from 
KPNO data. The tick marks indicate the phases of observation. Note that the longitude runs in the opposite sense to phase, 
with longitude 0° at phase 1.0 and longitude 360° at phase 0.0. Right-hand panel: The mean fractional spot occupancy of the 
reconstructed image as a function of latitude. 




Fig. 5. Same as Fig. [IJ but for the secondary component (G5 V) of the HR 1099 system. 



by using telluric features fsee lCollier Camerorll999h . This 
ensures reasonably accurate radial velocity determina- 
tions. Taking into account the different phase ranges and 
different instruments used in both data sets, one would 
expect some spurious scatter in the radial velocity. 

3.6. Final Image reconstruction 

Surface images of HR 1099 were obtained from the 
spectroscopic deconvolved line profiles and photomet- 
ric light curves using the DoTS surface imaging code 
ijCollier CamerorJl997l). Based on the M EMSYS algorithm 
developed bv ISkilling fc Brvanl l|l984lh DoTS is a maxi- 
mum entropy code which is used for the Doppler tomog- 



raphy of stellar surfaces. ICollier Cameronl ( 1997t) describes 
the stellar surface model and the integration scheme. 



Although some authors (e.g. iDonatil Il999(l have ne- 
glected the contribution of the G5V secondary compo- 
nent to the global activity of the system, we do try to 
obtain an image from both components of the binary sys- 
tem. The maximum entropy fits to the KPNO and MSO 
spectroscopic data sets are shown in Fig. We observe 
that the deconvolved profile of the system HR 1099, is 
well fitted at almost all the orbital phases. However, for 
phases near 0.0 and 0.5 the fit is not as good due to the 
degeneracy introduced by the overlapping profiles. Fig. |21 
shows the maximum entropy fit to the photometric light 
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Fig. 6. Maximum-entropy-regularized image reconstructions, for the primary component (Kl IV) of the HR 1099 system, from 
MSO data. The tick marks indicate the phases of observation. Note that the longitude runs in the opposite sense to phase, 
with longitude 0° at phase 1.0 and longitude 360° at phase 0.0. Right-hand panel: The mean fractional spot occupancy of the 
reconstructed image as a function of latitude. 
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Fig. 7. Same as Fig. |SJ but for the secondary component (G5 V) of the HR 1099 system. 



curve. Figs. 14161 shows a clear high-latitude spot between 
phases 0.55 and 0.85. 

The maximum entropy reconstructed images for both 
KPNO and MSO data sets, Figs. 14171 are a combination of 
the spectroscopic and photometric data sets. Both KPNO 
and MSO data sets have been treated independently. A 
stellar surface grid with 90 latitude bands was used for 
the image reconstructions. 

As we have already said in §3, the spectroscopic re- 
constructions have poor latitude discrimination at latitude 
<30°. On the other hand, the photometric data is sensitive 
to structures situated at low latitudes, although at much 
lower resolution. All this implies that we have to weight 
the spectroscopic and photometric data sets in order to 



combine them properly. The DoTS surface imaging code 
includes a weighting factor, (3, that determines which data 
set, spectroscopic or photometric, will contribute most to 
the image reconstruction. The weighing factor can be set 
from (100% photometric data) to 1 (100% spectroscopic 
data). However, (3 is larger than 0.9 in most of the cases, 
due to the higher amount of information contained in the 
spectroscopic data. In our case, the value for (3 was set to 
0.90 and 0.92 for the KPNO and MSO data sets, respec- 
tively. 
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4. Results and Discussion 

A maximum entropy reconstruction of the spot distribu- 
tion on both components of the binary system HR 1099 
is shown in Figs. 14171 The mean fractional spot area as a 
function of latitude is also plotted for each map. It shows 
that medium-high latitudes are the main regions of spot 
coverage for the primary component. The secondary com- 
ponent shows a medium-low latitude spot around orbital 
phase 0.7. However, the are also other features on the sec- 
ondary component which are not reliable as we will explain 
below. 

There seems to be a global similarity between the 
maximum entropy reconstructed images for the primary, 
(Figs. 01 and EJ), where starspots are mainly observed at 
latitudes higher than 30°. In both maps the surface ap- 
pears to be heavily spotted at the phases 0.6 to 0.9, 
while spotted at the phases 0.2 to 0.6. All these spots, 
detected from KPNO and MSO spectroscopic data sets, 
are in agreement with those detected from photomet- 
ric s^a^jnodejbn^josmgcontemporaneous photometry, 
bv lOarcia- Alvarez et al.l l)2003h . These authors, based on 
spectroscopic data sets obtained from several sites dur- 
ing the MUSICOS 1998 campaign, also observed rota- 
tional modulation and chromospheric flares occurring at 
the same phases. 

Although our sampling does not cover all the or- 
bital phases, the low inclination of the stellar rotation 
axis (40°), indicates that most spot features that ap- 
pear at high latitudes are in view all the time, con- 
tributing to every line profile, independently from the 
phase. Both Figs. 0] and recover a structured polar spot 
which is off-centered, with r espect to the stellar rota- 
tion pole. iDonati et al.l <|2003l) reported very similar spot 
distribution on the primary component, from observa- 
tions around the same epoch, reconstructed from Zeeman- 
Doppler imaging spectropolarim etry. This near-polar spot 
was also found by other authors JVoet et al 1 ll999HDonatl 
Il999t IStrassmeier fc Bartusl l200o|l bv using independent 
codes. lHatzes et alJ l|1996j) presented a detailed study of 
the reality of polar spots. Their modelling simulation ar- 
gues strongly against gravity darkening, differential rota- 
tion, limb brightening, equatorial bright bands, chromo- 
spheric line filling, and unknown effects in line radiative 
transfer physics as possible causes of the spurious polar 
spots. They suggested an inclination dependence of the 
flattening of the core of the line profiles in RS CVn's as 
an evidence for the reality of the large cool polar spots. 
The long vertical feature observed in Fig. 0] at phase 0.6 
may be real, because the sampling around this phase is 
large enough. Although the entropy criterion expresses 
the uncertainty in latitude by smearing the image north- 
south at low latitudes, one has to bear in mind that the 
combination of the spectroscopic and photometric data 
sets obtaining the maps should have partly removed this 
problem. No low-latitude features were obtained in the 
reconstructed images. 



There are, however differences between the maps, 
which may be partly due to the different spectral reso- 
lution of the KPNO and MSO spectrographs. The higher 
resolution of the KPNO data (R=65 000) allow us to de- 
tect isolated features and groups of features, while in the 
lower resolution MSO data set (R=35 000) those features 
will appear as a single big spot (e.g., features at latitude 
60° and phases 0.7-0.9). The resolution problem should 
not be affected by the input S/N, because it is similar for 
both KPNO and MSO data sets. Note that we are con- 
fident that the data has not been over-fitted. This addi- 
tional fact could add spurious surface features. We have to 
take into account that the reduction process is never per- 
fect, combined with incomplete phase coverage. However, 
by stopping the iterations before the noise is being fitted 
we ensure a reliable result. 

The images for the secondary components obtained 
from both KPNO and MSO data sets show a few medium- 
low latitudes features. We observe a low-latitude spot 
around orbital phase 0.7 in the images from both data 
sets. This spot seems to mirror the structure seen on the 
primary component. However, one has to bear in mind 
that most of the evidence for the structure on the G5 V 
component comes from changes in the EW of the star's 
contribution around the orbit, since its profile is so nar- 
row. We performed a series of numerical experiments on 
our HR 1099 data sets, altering the effective temperature 
of the spots, in order to see if the spot distribution on the 
secondary is sensitive to the spot temperature assumed. 
We conclude that the spot observed on the secondary com- 
ponent is not affected by changes in the surface brightness 
of the spots on the primary. This fact confirm the relia- 
bility of our results, insofar as sensitivity of the image 
to the spot temperature is concerned. However, this low 
latitude feature observed on the secondary component, is 
in disagreement wit h the convect i ye ove rshoot dynamo 
models reported by iGranzer et alJ l)2000h , which predict 
the emergence of magnetic field at higher latitudes. We 
also observe two features which only appear in the MSO 
data set. They occur at phases where no observations were 
made an d may be ghost features due to a poorly sampled 
data set (jStout-Batalha fc Vo g tllT99l. 

The mirroring effect, observed on the HR 1099 binary 
system, might suggest that tidal forces may influence the 
spot distribution. However, our observational data set only 
covers a few orbital rotations and, in order to claim this 
mirroring as a real physical effect, it would be better to 
confirm it by observing the HR 1099 system during con- 
secutive years or eve n try to find similar beh aviours in 
other binary systems. lHatzes fc Kursterl 1)19991) reported 
the surface images on both components of the active 
pre-main-sequence binary V824 Arae (HD 155555) ob- 
taining mirror images, as it happens in our observa- 
tions. Their images look nearly identic al if one r e verses 
or fl ips one with re s pect t o the other. IPiskunovl (^996) 
and IPiskunov et alJ l|200l|) suggested that the presence 
of spots at sub-stellar points, obtained through Doppler 
imaging, facing the other component on the binary system 
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ER Vul, could be attribute to the reflection effect. This is 
unlikely to occur in our case considering the separation of 
the two components of HR 1099, ~11 Rq. Even more, the 
maps obtained for the secondary component only show 
a small spot (consistent in both images) compared with 
the much larger structure obtained for the primary com- 
ponent. iGunn et al ] lll997f) presented extreme ultraviolet 
and radio observations of the RS CVn system CF Tucanae. 
They detected a modulation in the EUV and radio flux 
levels which led them to suggest the existen ce of an intra- 
binary region of activity. IGunn et a l. ( 1999) found a sim- 
ilar region on the Algol system V505 Sgr. These intra- 
binary regions would imply a joint magnetosphere and 
possible interaction between magnetic loops between the 
two components of the binary system. IGunn et alJ 1 19991) 
suggested that their observations could alternatively be 
interpreted as a surface feature on the later type star fac- 
ing the companion. 
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